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Abstract 



Classical XY spins on a two dimensional triangular lattice with antiferromag- 
netic interactions are reconsidered. We find that the Kosterlitz-Thouless tran- 
sition associated to the U(l) symmetry appears at a temperature 0.0020(2) 
below the Ising transition at 0.5122(1) associated to the Z2 symmetry. The 
Ising transition has critical exponents different from the standard ones. Using 
extensive Monte Carlo simulations for equilibrium and dynamical properties 
we show that the lack of universality observed in previous studies is due to fi- 
nite size corrections not taken account. Likewise the Kosterlitz-Thouless tran- 
sition has a critical exponent ij ~ 0.36 larger than the corresponding standard 
value 0.25. Also the helicity jump at the critical temperature is smaller than 
in the ferromagnetic case in disagreement with theoretical predictions. We try 
using the concept of an "quasi Universality class" to reconcile the standard 
critical behavior observable at higher temperatures with the different quasi 
universal one close to the critical region. 
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I. INTRODUCTION 



The improvement in micro-fabrication has increased greatly the experimental studies 
of the Josephson- Junction arrays of weakly coupled superconducting islands. [0,0 Phase 
transitions in these arrays are similar to those in two dimensional XY spin systems and 
the application of a magnetic field introduces additional frustration effects. This explains 
partly the revival of interest in studying 2D frustrated XY spin systems. The theoretical 
problem is connected to the presence of two symmetries and their coupling. Indeed the 
continuous XY symmetry leads to a Kosterlitz-Thouless transition driven by the unbinding 
of vortex-antivortex pairs, while the frustration introduce an additional Ising symmetry. 
Also for Helium-3 film an Ising symmetry due to the p-wave order parameter exists besides 
the continuous symmetry of a phase like for Helium-4 films. 0-§[ Contrary to the three 
dimensional case, the critical region in films could be accessible to experimental observations. 
The model we are going to analyze has also a link to the physics of early universe, J7| where 
more complicated couplings between symmetries are considered, 

Extensive research has been done on two dimensional frustrated spin systems like the 

17] or related Zig-Zag models, [0,0 the triangular model 
J 2 model, mm the XF-Ising model [g|,gl and the 



Fully Frustrated XY model 
we are discussing 
Villain model. 



17,20-22 the J x 



Also other systems have the same symmetries: the 19 vertex model 
the ID quantum spins, 29 the Coulomb gas representation, p0|-|33|| the XY 



XY model. 



5^ , |5o| or the RSOS model []3B| Therefore it is of great interest to compare the results for 



these different models in order to find out whether a universal critical behavior exists. 

A priori only two possibilities exist for the critical behavior. First, if the two phase 
transitions connected to the two symmetries are at the same temperature, then a new 
universal behavior should result. Secondly if the transitions are at two different temperatures 
one could expect that the transitions are of the standard Ising and Kosterlitz-Thouless 
type, especially if the transition temperatures are further apart. The presence of topological 
defects or vortices could change the reasoning. Indeed numerous studies show transitions at 
different critical temperatures, however with exponents for the Ising transition which vary 
from model to model and therefore cannot belong to the standard nor to one universality 
class. The reason could be the existence of a line of fixed points generating the various 
exponents found in numerical works. A second proposition assumes only one new fixed 
point with the variation of critical exponents due to improper finite sizes corrections. |28 



Thirdly Olson suggested that a large screening length around the critical temperature of the 
Kosterlitz-Thouless (KT) transition prevents to see the true standard Ising behavior. |27| 

In this article we will show that several studies which indicate a non-universality of the 
exponents are due to finite size corrections. We are therefore led to favor the existence of 
only one fixed point. Moreover to reconcile the result of studies done in the finite size scaling 
region very close to the critical temperature with those done at slightly higher temperatures 
we will discuss introduce the concept of an "quasi universality class" and give a physical 
interpretation as function of the size of vortex. 

The behavior of the Kosterlitz-Thouless (KT) transition have been less studied numer- 
ically. Especially the critical temperature is difficult to find even for simpler ferromagnet 
XY system. We propose a new way using the Binder parameter to overcome this difficulty 
and we obtain a critical temperature 0.002 below the Ising transition. We obtain a helicity 
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jump smaller than the jump of the ferromagnetic system and an exponent 77 m 0.36 that is 
larger than 77 = 0.25 for the ferromagnet, contrary to the theoretical predictions. [30- 33 



Also for the KT transition there is a discrepancy between simulations at or near the 
critical temperature and at high temperatures. Indeed in this last region the exponent 
T) ~ 0.22. We introduce a "quasi universality class" to explain this crossover. 

A combination of the Metropolis algorithm and over-relaxation algorithm reduce the 
CPU time by an order of magnitude than the metropolis alone. With longer simulations, 
up to 10 times compared to previous studies, we gain two order better statistics. 

For the determination of critical exponents we have used the finite size scaling method. 
However, since the possible presence of a screening length could prevent to see the "true" 
behavior, |27| we have also utilized the properties of the system in the short time critical 
dynamic which allows us to verify our results. Both methods agree. 

The outline of the article will be the following. Section II is devoted to the presentation 
of the model. The Ising and the Kosterlitz-Thouless transitions are studied in section III 
and IV respectively, discussion and conclusion are disclosed in the last section. To avoid 
repetition with a previous article [El] we will refer often to it. 



II. MODEL 

We study the XY spins on triangular lattices with antiferromagnetic interactions. The 
Hamiltonian is given by: 

H = J y £S i .S j (1) 

(ij) 

= jJ2cos(e l -e,) (2) 

m 

where S» is a two component classical vector of length unit, J is the antiferromagnetic 
coupling constant ( J > 0), 9 varies between and 2ti and (ij) are the next nearest-neighbors. 

The competition between the interactions gives the famous "120°" structure where the 
spins are not collinear (fig. |l]) and where the frustration is divided amongst all links. Due 
to this structure, the simulated lattice sizes must be a multiple of 3. We simulate the sizes 
L =12, 18, 24, 36, 48, 60, 81, 105, 123, 150. Two ground states exist, not related by a global 
rotation, and therefore, in addition of the symmetry U(l) from the continuous aspect of the 
spins, an Ising symmetry is present. 

To compute the order parameter for the U(l) symmetry we divide the lattice in three 
sublattices (s = 1, 2, 3) with only parallel spins in the ground state. After having calculated 
the magnetization of each sublattice M s we sum them to obtain M: 

M = ^E|M S |. (3) 

N = L 2 is the total number of the lattice sites. 

The order parameter k of the Ising symmetry is the sum of the chiralities K{ of each cell 

(fig-H)- 
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"i = ^ [S! x S 2 + S 2 x S? + x S, 1 ] , (4) 

K =™IE^I> ( 5 ) 

where the summation is over all cells and N' = 3N is their number. The chirality of one 
triangle is parallel to the Z-axis and equal to ±1 in the ground state only. We note that we 
could take as a definition the sum of the signs of chiralities. The result should be similar 
because the length of the chiralities is not relevant for the critical behavior. 

To compute the critical properties we have to define for each temperature the following 
quantities: 

Xf = W (6) 



iV((K 2 )-(.) 2 ) 

v- - N(k2) (8) 

1? = M _ {E ) (9) 

^ = {J w ~ {E} (10) 

v? = ^ - <*> (id 



^ 1 - (12) 
T -vf-7^ <l S sm( "--^ |2> (14) 

i±7 is the energy, x is the magnetic susceptibility per site, V\ t 2 are cumulants used to obtain the 
critical exponent v, U are the fourth order cumulants, = Xi—Xj where is the coordinate 
of the site i following one axe, T is the helicity [^lj] corresponding to the increment of the 
free energy for a long wavelength twist of the spin system, [0,^] (...) means the thermal 
average. 



III. ISING SYMMETRY 

This chapter is devoted to the Ising symmetry. The first part is related to the properties 
in equilibrium, i.e. we calculate the various quantities after a time tth to thermalize the 
system much greater than the correlation time r. The average is done on a time t av which 
is also much greater than r. In a second part we will use the short time critical dynamic 
recently introduced, i.e. the reaction of the system to a quench at the critical temperature 
from an initial state. No thermalization is used (tth = 0). 
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A. Equilibrium properties 



1. Algorithm 



As explained in our previous article, |2J] we use a combination of Nmet Metropolis steps 
and Nor over-relaxation steps. [39~] The over-relaxation algorithm is "microcanonical" in the 
sense that the energy does not change under a step. This algorithm reduces considerably 
the autocorrelation time. We have thus two parameters (Nmet and N r) to fit in order to 
minimize the CPU time fcpu- Our implementation for the over- relaxation algorithm is six 
times quicker than the Metropolis algorithm: 

fcpu (N m et i N r) = r (N MET + — ^) (15) 
= rN MET (l + -^p-) (16) 

IV MET 

In Fig. |2| we have plotted the autocorrelation time r for the chirality n, at the critical 
temperature T = 0.5122 calculated below, multiplied by N M et, as function of N r/N M et 
in a log-log plot for a lattice size 36. r is calculated with the method explained in the 
Appendix of our previous article. The data are well described by: 

AT ( No * ^ (MS 

TN M ET = a L [— (17) 



N 



MET 



with Gt36 = 160 and 636 = —0.57. 

Using (|i~7|) it is not difficult to show that the minimum of fcpu dH) occurs for: 



Nmet 1 + b 

which is 7.95 for L = 36. We have determined b for several lattice sizes L. If we divide 
the result by L, we obtain the ratio cl- c\i ~ 0.18, C24 ~ 0.28, C36 ~ 0.24, ~ 0.21, 
cqo ~ 0.21. This ratio is nearly constant whatever and L are. This is in accordance 

with a conjecture of Adler [fHJ stating that this ratio is proportional to the correlation 
length, i.e. the size in the finite size region where we have done the simulations. We have 
chosen Nmet = 1 and Nor ~ 2L/9 for the simulations. 

We show in Fig. |3| fcpu fo r the chirality as function of the size of the lattice L in a log- 
log plot for various algorithms at the critical temperature T = 0.5122. A similar behavior 
is obtained for M 2 , i.e. the order parameter associated to the U(l) symmetry. A is the 
metropolis algorithm alone, Ai is used in combination with one step of an over-relaxation 
algorithm, A c l in combination with 0.22L steps of an over-relaxation algorithm, while Aw 



correspond to the use of the Wolff cluster algorithm. |4l 
The slope of A gives us the dynamical exponent 

z K = 2.30(4). (19) 

We will discuss the value of this exponent at the end of this section. If we compare fcpu 
for A and A cE we observe that the gain is about 10, which means that for the same time 
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of simulation we obtain ten times better statistics (i.e. "independent" data). We note that 
Wolff's algorithm is less effective than the Metropolis algorithm. This is understandable 
because this algorithm uses only one link at each time to construct the cluster and we know 
that three links must be taken into account (at least one cell), therefore the algorithm can 
not generate the "good" cluster. Indeed each cluster has about 80% of the sites of the lattice, 
which is too high. Moreover even if we would be able to construct a good cluster, there is 
no guarantee that the method to flip the spins inside would be efficient. B2" 



2. Errors and details of the simulation 

We follow the procedure explained in the Appendix of our previous article. 0] 
We use in this work the histogram method which allows us from a simulation done at 
T to obtain thermodynamic quantities at T close to T . However to reduce the systematic 
errors we do not save histograms for (M) , (M 2 ) ... as function of energy but save the data, 
that is the energy E, the magnetization M and the chirality k. To avoid the use of a large 
space on the hard disk the data is saved only every r s sweeps (see table |). This method 
slightly increases our statistical errors (approximative of 20%). However the systematic 
errors decrease and we have a better control of the total errors. 

Since the data of a Monte-Carlo simulation are not independent, the calculation of errors 
must be done carefully. For simple quantities like the magnetization, the calculation is done 
with the standard formula of statistic but with the number of "independent" data ti n d 
equal to ^1°^ = ^ A _[ c . [Ell Nmc is the number of Monte Carlo steps. However for more 

ZT J T s -pi ZT-\-T s l| ,J 



complicated quantities we have to consider the correlations between the components of the 



formula and use, for example, the jackknife procedure. |45| Formula (A8-A14) of our previous 



article can be used but we have to change N MC to N mc /t s and r to t/t s . However the 



error for the helicity is not given. The helicity can be written: 

T = (A) - (B) (20) 



where A and B are given by ([14]). Applying the same method as in Refs. [24], it is not 
difficult to obtain: 



AT 2 = ^^[(A 2 ) - (A) 2 + <B 2 > - (Bf - {(AB) - (A)(B))\ . (21) 

The simulations have been done using sizes between 12 to 150. In the table [I] we gave 
some details of the simulations where t t h is the number of Monte Carlo steps (i.e. one 
Metropolis step followed by 2L/9 over- relaxation steps) to thermalize the system; t av is the 
number of Monte Carlo steps to average. The third column gives the autocorellation time, 
the fourth the time between two consecutive measures r s . The last column is the number of 
"independent" data. We report the last line in table [TI|, it should be compared to previous 
studies. Our statistic is two orders greater than previous studies for similar sizes. This 
allows us to obtain better precisions for the quantities, typically one or two orders smaller, 
and therefore the critical exponents are more reliable. In particular we will see that the 
finite size corrections are important and this explains the variation of exponents found in 
previous studies (see discussion at the end of this section). 
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3. Results 



Our first task is to find the critical temperature T C K . The most effective way is to use 
Binder's cumulant (|T3| ) in the finite size scaling (FSS) region. We record the variation of 
U K with T for various system sizes in Fig. |] and then locate T C K at the intersection of these 
|jj since the ratio of U K for two different lattice sizes L and V = bL should be 1 



at T C K . Due to the presence of residual corrections to finite size scaling, one has actually to 
extrapolate the results taking the limit (ln6) _1 — > in the upper part of Fig. |5|. We observe 
a strong correction for the small sizes. However for the biggest sizes the fit seems good 
enough and we can extrapolate T C K as 

T C K = 0.5122(1) . (22) 

We note from the figure that it is a lower bound. The estimate for the universal quantity 
U£ at the critical temperature is 

U? = 0.632(2). (23) 



This value is far away from the two dimensional Ising value JJ^ sm9 ~ 0.611 |^7[ which is a 
strong indication that the Universality class associated to the chirality order parameter is 
not of Ising type. We note that it is not compatible with the value of the J\ — J2 model 
0.6269(7) |24| which, in the standard formulation, means that the two systems belong to 
two different Universality classes. However since this transition is a coupling between two 
it is not certain that this quantity stays universal. 

At T = T C K the critical exponents can be determined by log-log fits. We obtain v K 
from V* and V 2 K (Fig. |), 7 K jv K from \ K an d X2 (Fig- 0)> an d (3 K /v K from k (see Fig. ||). 
We observe in these figures a strong correction to a direct power law. It is worth noticing 
however that X% shows smaller corrections. Using only the three (four for xT) largest terms 
(L =105, 123, 150) we obtain: 

v R = 0.815(20) (24) 
7 7z/ K = 1.773(9) (25) 
p*/u K = 0.110(6) . (26) 

The uncertainty of T C K is included in the estimation of the errors. The large values of 
our errors are due to the use of only few sizes for our fits. The values obtained for four 
consecutive sizes are interesting. For the exponent v we obtain: ^12-18-24-36 — 0.909, 

^18-24-36-48 = 0.909, ^24-36-48-60 = 0.903, ^36-48-60-81 = 0.884, ^48-60-81-105 = 0.862, 

^60-81-105-123 = 0.835, ^8i-i05-i23-i50 = 0.822. These values cover a large range of the 
data obtained by previous studies (see table ||) and we strongly suspect that the lack 
of Universality (at least in the critical exponents) is related to the corrections previously 
ignored. In the conclusion of this section we will show which important physical informations 
we are able to obtain from the sign of these corrections. 

We have tried to introduce a correction to calculate the exponents, for example for V* = 
(1 + L _wK )L 1 / !/K , we obtain oj k = 1.2(5) and values for critical exponents fully compatible 
with 
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The values given in (|24| - p6D use the properties of the free energy at the critical temper- 
ature. But an error on T C K leading to an error on the exponents, it is therefore advisable 
to find them without the help of T C K . This can be done using the whole finite size scaling 
region with the method given in. jP| It consists to plot, for example, the susceptibility 
function of U K , choosing the exponents so that the curves collapse. This fit is 
more reliable than the fit at the critical temperature since as it does not depend only on the 
results at T C K but on a large region of temperature. However the errors are a little bit larger. 
We show in Fig. |9|-|TT| the results for three choices of 7 K /V K 1.82, 1.78 and 1.74. Obviously 
the second is the best and we obtain 7 K jv K = 1.78(2) compatible with the result of (P5p. 

In Fig. [12] we have plotted the exponent rj K = 2 — 7 K /V K calculated with a direct fit of the 
susceptibility (like in Fig. |7|) but changing the temperature. In the same diagram we have 
plotted the result (A77") found by the fit of \ K as function of U K . This gives an estimate of 
the critical temperature (see figure) compatible with T C K . 

The results are summarized in table |T|. In the next section we will try to compute the 
properties of this transition using the short time critical behavior. This method does not 
suffer from the same problems as the finite size scaling method and allows us to check our 
results. Finally we will compare our results with those of previous studies. 



B. Short time critical behavior 

In this section we want to study our system using the dynamical behavior near the 
critical temperature. This particular method consists of quenching the system from zero 
temperature to the critical temperature and studying the short time dynamics, i.e. before 
reaching the equilibrium properties. It seems strange that the system shows critical behavior 
although it has not reached equilibrium. Indeed for a long time this fact has not been utilized 
in numerical simulations. However it has been demonstrated that between a time t{ and a 
longer time tf a region exists where the calculation of the critical properties of the transition 
is possible [49-Q (see the review of Zheng ||52|| ). In this region the system has lost the non- 



universal informations (t > tj) but the correlation length £(t) is much smaller than the 
correlation length £(oo) of the system in the equilibrium phase. Since we are at T c the 
correlation length £(oo) is infinite but in our simulation it is upper-bounded by the size L 
of the lattice and the condition written as £(t) -C L. After tf this condition is not satisfied 
and the system shows a crossover to equilibrium properties. In the region tj < t < tf the 
system shows very simple power law even for the Kosterlitz-Thouless transition. 



The quantities to compute are the same as in formula fl3|-|l"3|) but the average (. . .) is now 
taken on different realizations, i.e. starting from the ground state with different random 
numbers for the Metropolis procedure. 

To determine the range of time and the size L that we use in the simulations we have to 
find tf where the system begins to show a crossover to the equilibrium properties. In Fig. |13 



we show (ft(t)) as function of the time t for different sizes L at the critical temperature 
T = T C K = 0.5122 in a log-log plot. We observe that for small systems t/ ~ 100. It grows 
for larger systems and since the sizes L = 201 and L = 300 give similar results for (n(t)) in 
Fig. O tf is bigger than 10 4 . Therefore we can use the entire region to calculate the critical 



exponents for the greatest size. 
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We have studied systems with size L = 300 for a simulation time t up to 10,000 and 
we have averaged over 6000 realizations. This reflects more than one order statistics better 
than previous studies on the fully frustrated XY model (FFXY). For the algorithm we 
use only the Metropolis algorithm (Model A in the classification of Ref. f54|). Here we do 
not use the over-relaxation algorithm. The errors are calculated dividing the data in three 
sets of 2000. 

We want now to verify the critical temperature T C K = 0.5122 found in the previous 
section. For this purpose we have simulated our system for three different temperatures at 
and around this value. We know that at the critical temperature the magnetization shows a 
power law behavior, a straight line in a log- log plot. For T^T C the system shows important 



corrections as observed in Fig. |14 . 



With our knowledge of the critical temperature found above, we are able to calculate 
very precise exponents. In the short time critical dynamic the quantities (|3|-|T3D have at T c 
the behavior: 

(«(*)) oc r^/*""*") (27) 
X K oc (28) 
V k K t i/(r"*") (29) 

U R oc t d/zK (30) 
where z K is the dynamical exponent and d the dimension of the space, i.e. d = 2. 



We first compute z using (|30|) . In Fig. |15| we have plotted U K as function of t. We see 



that from £j ~ 100 the curve shows a linear behavior. Using a fit from t = 300 to t — 10, 000 
(shown by the two arrows) we are able to obtain 

z K = 2.39(5). (31) 

This value is consistent with the value found in equilibrium properties (z* = 2.30(4)), but we 
are more confident in the result of dynamical properties which are less subject to systematic 
errors (see Appendix of Ref. ||24| ). It is also consistent with the value of z — 2.31 — 2.36 found 



studying the equilibrium properties of the Fully Frustrated XY model. |55| Nevertheless it 
is in contradiction with the value z K = 2.17(4) found studying the dynamical properties of 



the last model. |L6| However a slight change of the critical temperature could have a strong 
influence on z and the calculation have been done with more than one order less statistics 
than our work. 

To obtain the other exponents we could use the formula (^-^), however we prefer to 
obtain directly v K and 7 K plotting U K and x K as function of V* and r] K = 2 — / ~f K / v K plotting 
X K as function of U K . This has been done in Fig. [16] and Fig. |17]. Fits have been done for 
300 < t < 10, 000 represented by the arrows. We have also calculated /3 K plotting (k) as 
function of V" (not shown). We obtain: 

v K = 0.818(9) (32) 

7 K = 1.445(20) (33) 

(3 K = 0.0967(13) (34) 

if = 0.235(7) (35) 
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Our results are summarized in table |T|. The exponents calculated by the two methods 
(equilibrium and dynamical) agree very well. However we would prefer the results from 
short time critical dynamics since they are less sensitive to the finite sizes corrections. 



C. Comparison with previous studies 

Before discussing the interpretation of our results, we will compare them to those found 
previously on various models. 

In table |lj we review the numerical studies related to our work. We think that the errors 
quoted in this table are mostly too optimistic estimates. 

We first compare our work with previous simulations using the same methods as this 
work (Monte Carlo Finite Size Scaling-MC FSS). There are two important columns: the 
maximum size L max used (fourth column) and the number of independent data ti n d (seventh 
column). Indeed we have seen that strong finite size corrections exist in our system and we 
cannot see them if the largest size is too small. When comparing the calculations with a 
similar largest size (L ~ 150), |H||2"2"| , |2"5[ we observe that our statistics (t ini i) is between 100 



to 1000 times better. This allows us to see the corrections not seen before and to explain the 
lack of Universality observed in previous studies (see subsection "equilibrium properties"). 
The results, for the triangular lattice, in favor of the Ising ferromagnetic Universality class 
(u = 1) f22| are not reliable due to too small statistics {t in d = 15). 

For the other methods none of the authors have reported the strong corrections which 
should exist. Nevertheless almost all results are compatible with ours iy ~ 0.795 — 0.82) 
and are in favor of a single Universality class. Only results on the XY-Ising model seem 
to show different Universality classes but simulation sizes are so small that results are not 
conclusive. 

Our results are thus strongly in favor of the picture of p8| which suggests the phase 



diagram drawn in Fig. [18|. A and C are free parameters. Similar phase diagrams appear for 



the Ji — J2 model, fiZM the Zig-Zag model fli8|,|T9|1 and the RSOS model [3E|. A variation in 



free parameters changes only the initial point of the renormalization group flow on the line 
PT but all trajectories converge to the same fixed point F, i.e. a single Universality class. 
The finite size corrections will be more or less important following the initial point. To verify 
this interpretation it would be useful to test several initial points of the line (PT) of the 
XY-Ising model with similar sizes as in this work (L max ~ 150) to see if the systems belong 
to the same Universality class and, maybe more interesting, to observe the corrections to 
scaling. The picture should be very similar to the Potts model with disorder where this 
crossover, and therefore the corrections, have been observed. KPI] One other numerical study 



also interesting should be the calculation of the properties of the 19-vertex model by the 
Density-Matrix Renormalization Group (DMRG). The model has already been studied by 
Monte Carlo Transfer Matrix PS| but only for small sizes (15 2 ). On the contrary the DMRG 
allows to treat L chains of infinite sizes. We note that the method has been proved to be 
very efficient for the 19-vertex for the ferromagnetic system (i.e. with different internal 
parameter). |K| 

We now discuss the results of the Monte Carlo (MC) simulation in the high temperature 
(HT) phase. In this region we have to keep the correlation length £ much smaller than the 
size L of the lattice: (<L. studies on the FFXY model have been done by Nicolaides ITT 
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and by Jose and Ramirez. [12| They found two very different results. However with a close 
look on their simulations we have found that the results of Nicolaides were not reliable and 
his errors more than strongly underestimated. Therefore the HT MC seems in favor of a new 
Universality class but with an exponent v ~ 0.89 larger than found in the finite size scaling 
MC {y ~ 0.81). On the other hand Olson |27j has done a HT MC on the Villain model, 
i.e. where the spinwave are decoupled from the vortex, and has observed a very interesting 
behavior. He proposed that a screening length A = C,kt exists, due to the Kosterlitz-Thouless 
transition and, that in addition to the condition C,i S i ng <C L, the sizes have to satisfy the 
condition L. In this case the system shows a standard Ising behavior v = 1 while it 

is in the crossover to this behavior if the condition A <C L is not respected. Indeed he found 
some indication that the exponents v ~ 1 when the two conditions are satisfied. However 
since the numerical results are not extremely precise, further simulations are needed to prove 
this interpretation. 

Let's assume that this last interpretation is correct. It does not prove that there is no 
new fixed points contrary to the claim of Olson. Indeed we have shown in this work and 
in |SJ that the corrections lower the result for v (from 0.90 to 0.80) and therefore that our 
system is not in the crossover to the ferromagnetic Ising behavior v — 1. We then conclude 
that a new "stable" fixed point exists for a range of size L. We call this new fixed point an 
"quasi fixed point" by similarity with the 3D case where an "almost second order" exists for 
a certain range of size L before showing the "true" first order transition. p6| , p7| This can 
be understood if we admit that the Kosterlitz-Thouless transition coupled the spins at large 
distance by the intermediate of £kt- Therefore it has tendency to bring the critical behavior 
of the Ising symmetry to the mean-field solution [y = 0.5) exact for infinite interactions. 
Since ^kt is finite above the critical temperature the thermodynamic limit of this "quasi 
fixed point" is not stable. 

To be exhaustive we mention the study of the transition W-Potts(g). |58| We have 
shown that the transition is of first order whenever q > 3. In this case the transition will 
stay of first order even in the thermodynamic limit, since the correlation length is finite. 
Therefore the new fixed point (or his absence) is stable. This fact is rather against the 
interpretation of Olson but it is possible that the stability of the new fixed point changes 
between q = 2 (Ising) and q = 3. On the contrary it would say that for the W-Ising model, 
the two behaviors in the finite size scaling region [y « 0.80) and far away of the critical 
temperature [y = 1) are different and stable in the thermodynamic limit. 



IV. U(l) SYMMETRY 

We will present in this section our results for the Kosterlitz-Thouless (KT) transition 
|61,62] associated to the U(l) symmetry. This transition is characterized by the unbounding 



of the vortex-antivortex pair at the critical temperature. At low temperature there are 
only some pairs of vortex-antivortex while at high temperature alone vortex (or antivortex) 
can exist. This picture is considered as correct in the ferromagnetic and frustrated case. 
However the two cases could show some differences for the jump of the helicity at the 
critical temperature T C M and for the exponent rj. 

The helicity T is the answer of the system to a twist in one direction. It has been shown 
that at T C M this quantity jump from to an universal quantity \E2\ for all ferromagnetic 
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systems. Moreover the behavior of T as function of the lattice size has been predicted. |63[ 



The use of T has been proved to be very useful for the ferromagnetic case and in particular to 
determine the critical temperature. |64| The situation is not so clear for frustrated systems. 
It has been suggested that the jump at T C M could be "non universal" i.e. different from the 
ferromagnetic case p0|-p5| and no scaling with the system size has been proved. Therefore 
rather than using the helicity to obtain informations on the transition, we will use a new 
method for the KT transition that we have introduced in. It consists in using Binder's 

cumulant to study this transition. It was proved in these articles that, contrary to the 
common belief, the Binder cumulant for ferromagnetic XY systems crosses for different 
sizes, allowing thus an estimate of the critical temperature and moreover an estimate of the 
exponent 77 without the precise knowledge of the critical temperature. The ferromagnetic rj 
has been proved equal to 1/4 p2l with logarithm corrections |J7| while it has been predicted 
to be smaller in frustrated cases. f30H35 



In this section we will first present our results for the equilibrium properties and second 
verify them using the short time dynamical properties. 



A. Equilibrium properties 

1. Critical temperature and critical exponent rj 

We proceed in a similar way as for the Ising order parameter, replacing k by M. We 
record the variation of U (|T2"D with the temperature for various system sizes in Fig. [19. 



We want to underline the differences between our result in the frustrated case and in the 
ferromagnetic case: the crossing region is one order of magnitude less than for the standard 
XY model (compare with Fig. 1 of Ref. f65|). We then locate the intersection of these curves 
and plot the results in the lower part of Fig. |5|. 

Let us first consider a power law behavior at T > Tc for this system. We then have to 
consider a linear fit for (lnfr) -1 — > 0. We observe corrections for the smallest sizes L =12, 18 
and 24 but the others seem to converge to the temperature 

T C M = 0.5102(1) . (36) 

We note from the figure that it is an upper bound for the critical temperature, i.e. it cannot 
reach T C K = 0.5122(1) and therefore we obtain two transitions. 

Secondly we consider the behavior to be exponential as in the standard XY model. In 



this case Fig. 2 of Ref. |65j shows that a linear fit could be wrong and that a " crossover" to 
a different critical temperature could be observed for bigger b, i.e. greater sizes. However 
contrary to the ferromagnetic XY model, the region of crossing is very small and the different 
linear fits tend only to one critical temperature. We observe the same behavior for the J\ — J2 



model. [23J We conclude that the linear fit works well enough and we will show below strong 
arguments in favor of the temperature 

With the help of the critical temperature we have found an estimate of U M at T C M fitting 
the value with a law U M = + aL-°: 

Uf = 0.6497(12) . (37) 
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The exponent r\ is obtained by a log-log fit of x¥ as function of L shown in Fig. 0. We 
obtain 

2 - rj M = 1.635(10) (38) 
V M = 0.365(10) (39) 

The fit has been done using only the four largest sizes (L = 81, 105, 123 and 150), disre- 
garding the smallest sizes which show small corrections. We note that for small sizes 7] M is 
smaller, for example if we use only the sizes from 12 to 60 we obtain 7] M = 0.33. Therefore 
if we take into account the corrections, the exponent moves away from the standard value 
0.25 for the ferromagnetic systems. This shows that we are not in a crossover to this last 



behavior. Our value is in agreement with those found for the J\ — J2 model. |24| Our results 
are summarized in table |TJ. 

From a theoretical point of view the KT transition has an exponential behavior, i.e. a 
correlation length of the form £ ~ exp[2? (T — T c )~ u ], however a power law behavior like 
(£ ~ (T — T c )~ u ) can not be excluded numerically. In the latter case the critical exponent 
v can be calculated with the cumulant V 2 M (|llD- We have obtained v M = 1.18(10) and 
the results show important corrections. This value being in contradiction with the value 
v = 0.92(3) found in the J\ — J2 model, we can exclude a power law behavior. 

As for the Ising order parameter, the calculation of the exponents has been done at 
the critical temperature but an error on T C M leads to errors on the exponents, it is then 
interesting to find them without the help of T C M . This can be done using the same method 



as described before. We have shown in Ref. |65| that this method is accurate enough in order 



to obtain 77 whatever the type of the behavior is (power law or exponential). In Fig. p0|-^2 
we show our results for three values of r\ M , 0.34, 0.375, 0.41. Obviously the second value is 
the best and obtain: 

r) M = 0.375(20) (40) 

which is compatible with (|39|). 

In Fig. |12| we have plotted the exponent i] M calculated with a direct fit of the susceptibility 
(like in Fig. [7]) but changing the temperature. In the same diagram we have plotted the 
result (Arj M ) found by the fit of x¥ as function of U M . This gives an estimate of the critical 
temperature (see figure) compatible with T C M . 

We now want to demonstrate that the two transitions cannot appear at the same tem- 
perature. First we have shown in Fig. |5]that T C K = 0.5122(1) is a lower bound for the critical 
temperature while T C M = 0.5102(1) is a upper bound. Moreover if the Ising transition ap- 
pears at T C M = 0.5102, the exponent rf must be equal to 0.115 (7" = 1.885) by a direct 



log- log fit at this temperature (see Fig. O). In Fig. 23 we have plotted \ K as function of 



U K for this value. The curves do not collapse. Now if the Kosterlitz-Thouless transition 



appears at T 6 K = 0.5122 then r) M = 0.461 (see Fig. [12]) and the curves X2 as function of 



U M do not collapse in Fig. ^4[ Therefore, even if the critical temperatures are different by 
only 0.4%, we are able to conclude that the KT transition appears at lower temperature 
than the Ising transition. Moreover we want to stress the fact that the phase diagram of the 



FFXY pi is not in contradiction with our picture. In this work the authors have shown 



that with varying a parameter (J'/ J) it is possible to decouple the two symmetries with the 
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Ising transition appearing at lower temperature than the KT transition. However in this 
case the Ising transition is due to the contraction to (or n) of the turn angle between the 
spins and not to the mixed of different chirality signs. |66| Therefore two physical properties 
exist near J' / J = 1 which lead to a complicated phase diagram with a cross of the critical 
lines. 



2. Helicity 



We present know our results for the helicity T (|T4]). In Fig. ^ we have plotted T as 
function of the temperature. T s is the temperature of simulation, T C M is our estimate for the 
critical temperature from above. 2T/ir is the universal jump for a ferromagnetic system. 

Tq is the best fit of T with the scaling form as function of the system size L valid in 



ferromagnetic systems |63 



21nL + c y 



(41) 



where Tj um p is the jump at T c and c a free parameter. We note that our results are very 
similar to those of Lee and Lee. 
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Since it is difficult to obtain, with the histogram 
method, informations for the biggest sizes too far away from the temperature of simulation 
(T s = 0.511) we take their result for Tq: 



T n = 0.501(1' 



(42) 



This value does not agree with T C M = 0.5102(1) found above. If we admit To as critical 
temperature the exponent r\ M must be equal to 0.22 (see Fig. [12]) and the curves x¥ as 
function of U M do not collapse in Fig. ^7|. Therefore we can doubt the validity of (41) for 
the frustrated case. 

More interesting is the comparison of T at T C M = 0.5102 with the ferromagnetic jump 
Tjump = 2T/7T. In Fig. |26| we have plotted these two quantities as functions of the size 
L. Our results suggest a jump at the critical temperature smaller than the ferromagnetic 
jump in contradiction of the suggestion of. p0|-[35|] This surprising result could induce some 
doubts about our previous results. In order to verify them, we will present now the dynamical 
properties of the KT transition. 



B. Short time critical behavior 

This analysis is very similar to those done for the Ising symmetry. Zheng and co-workers 



52| , |53H have proved that a system quenched from a zero temperature to the critical tem- 
perature T c shows that a KT transition has a similar picture as the Ising transition, i.e. a 
power law behavior as function of the time t: 

(M(t)) oc t-" M /( 22 M) (43) 

U M oc t d ' zM (44) 

If the system is not quenched to the critical temperature corrections are present and the 
behavior ceases to be linear in a log-log plot. 
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In Fig. ^ we have plotted (M) as function of t for different temperatures. We observe 
a linear behavior for T C M = 0.5102 and clearly a deviation from this behavior for the other 
temperatures and particular for T = 0.5010. This is a strong indication that T C M is the 
good choice of the critical temperature. 

We have calculated the exponents rj M and z at T = T C M from fl4~3| - f44"D by similar methods 
used for the Ising symmetry. Since the average has been done only on one sample of 500 
configurations we are not able to calculate the errors. We obtain: 

z M = 2.10 (45) 
r] M = 0.36 (46) 



The value of r] M is in agreement with those found in equilibrium (see table pl|) . 

In conclusion we found that, numerically, the Kosterlitz-Thouless in the Finite Size 
Scaling (FSS) region has a complete different behavior than for ferromagnetic systems. We 
found an exponent r] M ~ 0.36 larger and a jump of the helicity smaller in frustrated case. 

Our exponent i] M is in agreement with our results for the J\ — J2 model ^jf_j 2 ~ 0.35 
(see table |HI|) with a similar method as this work, and with the result found in the YYXY , 



Vffxy^o.34, using a transfer matrix method. [£q] However they are in disagreement with the 
values found in the high temperature (HT) region 77 = 0.20 — 0.25. |T2|j27| Although new 



and surprising, we believe that our results are reliable and present hereafter the numerical 
arguments in favor of them: 

1. Fig. [| is very suggestive for the value of the critical temperature and shows that the 
behavior of the Kosterlitz-Thouless transition in frustrated systems is very different 



from the ferromagnetic case. [24 



2. The behavior of \¥ as function of U (Fig. PU|-P3) which give rj 



-,M 



3. The Universality of our results for 77 for different models (see table [I I 



4. The dynamical properties shown in Fig. ^8] are an another indication that our choice 
for critical temperature is correct. Moreover the exponent 77 is in good agreement 
with those found in equilibrium properties. 

Since the two calculations (in the HT and FSS regions) seem correct we suggest an 
hypothesis similar to those given for the Ising symmetry, i.e. different behavior in the FSS 
and HT regions. We give hereafter a physical interpretation of this behavior. At high 
temperature the radius of the vortex are very small and they are slightly connected together 
while the correlation length for the Ising symmetry is bigger than the radius of the vortex. 
The situation is different in the FFS region. In this case the radius of the vortex are large 
enough and the properties should depend on the domain walls due to the Ising symmetry 
with the Ising correlation length smaller or of the same size as the radius of the vortex. 
This argument is exactly the opposite of those given to explain the difference of behavior for 
the Ising transition: the "true" ferromagnetic Ising behavior appears only when is much 
larger than a screening length A, i.e. the £,kt- This difference of argumentation reflects the 
differences between the two symmetries: the Kosterlitz-Thouless transition is a transition 
driven by topological defects (vortex-antivortex), i.e. by "local" behavior, contrary to the 
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Ising transition. We call our picture " quasi Universality class" in a similar way as the Ising 
symmetry. The two interpretations diverge in the sense that in the thermodynamic limit 
(infinite size) the new Universality class for the Ising symmetry could be not stable (however 
see Ref. ||58|| ) while it is for the Kosterlitz-Thouless transition. However the two behaviors 
are unstable in the high temperature region, i.e. sufficiently far away from the transitions. 



V. CONCLUSION 

This article is devoted to the study of the triangular frustrated XY spin system in two 
dimensions by extensive Monte Carlo simulations. 

This system is characterized by two symmetries: a U(l) symmetry due to the continuous 
nature of the XY spins and an Ising symmetry due to the degeneracy of the ground state. 
Our study is restricted to the region where the system sizes that we simulate are smaller than 
the values of the correlation lengths of the U(l) symmetry like those the Ising symmetry 
should have in an infinite system. 

We have shown that the Kosterlitz-Thouless transition connected to the U(l) symmetry 
appears at lower temperature than the Ising symmetry. Our very large statistics (two order 
larger than previous studies) allow us to take into account the corrections to the scaling 
laws. With the knowledge of corrections we have shown that the Ising transition belongs to 
a new Universality class, i.e. the renormalization group flow tends toward a new "stable" 
fixed point. These corrections explain also the lack of Universality found by several authors. 
Knowing that the transition has a different behavior in the high temperature phase we have 
introduced the idea of an "quasi fixed point" or "quasi Universality class". In this inter- 
pretation this new fixed point could be unstable in the thermodynamic limit but the true 
behavior cannot be reached in our limited accessible sizes in numerical simulations. Accept- 
ing this interpretation, the situation is similar to those who appear in frustrated systems in 
three dimensions. |)7j We note that our interpretation goes further than the interpretation 



of Olson [27] who predicts just an impossibility to obtain the true Ising behavior in the finite 
size scaling region. 

We have also studied the Kosterlitz-Thouless transition associated to the U(l) symmetry. 
We have shown that this transition belongs to a new Universality class characterized by an 
universal exponent r] ~ 0.36 and a helicity jump smaller than those in the ferromagnetic case 
in disagreement with theoretical predictions. We have introduced the concept of an "quasi 
Universality class" for this transition to explain the different behaviors observed between 
the finite size scaling and the high temperature regions. In this interpretation the new 
behavior is stable in the thermodynamic limit near the critical temperature but not when 
the temperature is much larger than the critical temperature. We have given a physical 
interpretation of this new concept. 

To bring new informations and verify our predictions for these systems we have three 
possible ways. The first way is to improve the theoretical approaches of these systems. In 
particular it should be extremely interesting to have a theoretical result for exponents of the 
new fixed point for the Ising symmetry. However since it is difficult to describe the coupling 
between the two symmetries, the precise calculation of exponents seems, for the moment, an 
impossible task. Another possibility is the experimental approach in the Josephson-junction 
array |T]H] or in films of 3 He. PH6J However one of the problems of experiments is to have 
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access to the interesting quantities. For example it is usually difficult to measure the chirality 
in these systems because it is not related directly to measurable quantities. The easiest way 
to obtain reliable information seem to be numerical simulations. Indeed the use of powerful 
computers and fast algorithms allows to obtain precise results out of reach some years ago. 
We have to study different models to observe the Universality class in the finite size scaling 
region but also in the high temperature region to verify the picture predicted in this work. 
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TABLES 



Uh 

3.10 5 
3.10 5 
5.10 5 
5.10 5 
5.10 5 
7.10 5 
7.10 5 
1.10 6 
1.10 6 
1.10 6 



30.10 6 
30.10 6 
30.10 6 
20.10 6 
18. 10 6 
30.10 6 
30.10 6 
27.10 6 
21. 10 6 
17.10 6 



11.9(1) 
13.0(6) 

50(2) 
73.(3) 
98(5) 
163(13) 
236(18) 
300(28) 
404(66) 



100 

"Too 
Too 
"Too 

100 

"Too 

~L63 
"236 
"300 
"404 



2t k .+t.s 

2.4 10 5 
2.4 10 5 
2.0 10 5 
1.0 10 5 
7.3 10 4 
1.0 10 5 
6.1 10 4 
3.8 10 4 

2.3 10 4 

1.4 10 4 



TABLE I. 
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TABLE II. W estimated using Ref. ||; W estimated using this work; ^^estimated using Ref. 



symmetry 



n 







Z 2 (eq) 



0.5122(1) 



0.6320(20) 



0.815(20) 



1.450(42) 



0.227(9) a 



0.0896(71) 



Z 2 (dyn) 



0.818(9) 



1.445(20) 



0.235(7)" 



0.0967(13) 



(eq) 



0.56465(8) 



0.6269(7) 



0.795(20) 



1.391(43) 



0.250(10) a 



0.101(10) 



U(l) (eq) 



0.5102(1) 



0.6497(12) 



0.365(10) 



(7(1) (dyn) 



[/ ( 1 )Ji-J 2 ( eq ) 



0.36 



0.56271(5) 



0.638(5) 



0.345(5) 



TABLE III. Calculated using 2 - rj = 7/zA 
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TABLE CAPTIONS 



Table [p Number of Monte Carlo steps to thermalize tth and to average t av as function of 
the lattice size. t k are calculated with shorter MC runs. We collect the data every t s . The 
last column gives the number of "independent" measurements. 

Table fTI|: Results for the two dimensional frustrated XY systems with a breakdown 
of symmetry U(l) ® Z 2 . Critical exponents are associated to the Ising symmetry. 1= 
Ising; FFXY= Fully Frustrated XY; CG= Coulomb gas; TA= Triangular Antiferro- 
magnetic; MC= Monte Carlo; HT= High Temperatures; FSS= Finite Size Scaling; Mi- 
cro=Microcanonic TM= Transfer Matrix; dyn= dynamic; RSRG= Real Space Renormal- 
ization Group; RSOS= Restricted Solid on Solid. The indices 1,2,3 for the Zig-Zag or the 
XY-lsmg model refer to different choice of internal free parameter. ti n d = NMc/^max] 

Table ITT: Summary of our results for the Ising symmetry (Z2) and the XY symmetry 
U(l). The results for the J\ — J 2 model is from Ref. [p4fl - (eq)=equilibrium properties; 
(dyn)=short time dynamical properties. 
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FIGURE CAPTIONS 



Fig. |l|: Ground state configuration for the triangular model. The chirality of each triangle 
is indicated by + or — . 

Fig. 0: Autocorrelation time r times the number of Metropolis steps N MET as function 
of N OR /N MET . N r is the number of over- relaxation steps. 

Fig. |[ CPU time (proportional to the autocorrelation time r) for the standard Metropo- 
lis algorithm (v4 -square) , in combination with one over-relaxation step (Ax-diamond), in 
combination with cL ~ 0.22 L over- relaxation steps (A C £-circle), and for the Wolff's cluster 
algorithm (A^-triangle) . 

Fig. ^: Binder's parameter U K for the Ising order parameter function of the temperature 
when L = 12 to 150. The arrow shows the temperature of simulation T s = 0.511. 

Fig. [|: Crossing T plotted vs inverse logarithm of the scale factor b = L'/L. The upper 
part of the figure corresponds to U K while the lower part to U M . We obtain T C K = 0.5122(1) 
and T C M = 0.5102(1) with a linear fit (see text for comments). 

Fig. H Values of V* and U 2 K as function of L in log-log scale at T C K . The value of the 
slopes gives \jv K . Strong corrections for small sizes are visible. Only the three largest sizes 
are used for the fits. The estimated statistical errors are smaller than the symbols. 

Fig. 0: Values of x K an d X2 as function of L in log-log scale at T C K and x¥ a t T^ 1 . The 
slopes give j/u — 2 — 77. Strong corrections for the small sizes for \ K are visible. Only the 
three largest sizes are used for the fit for \ K and the four largest for the fits for X2 an d X2 ■ 
The estimated statistical errors are smaller than the symbols. 

Fig. ||: Values of (k) as function of L in log-log scale at T C K . The slope gives /3 K /u K . 
Strong corrections for the small sizes are visible. Only the three largest sizes are used for 
the fits. The estimated statistical errors are smaller than the symbol. 

Fig. g x K L~ lK l vK as function of U R with 7" jv K = 1.82 for sizes from L = 24 to 150. The 
curves do not collapse in one curve. 

Fig. [TO]: x K L^ lK ^" as function of U K with 7 K /z/ K = 1.78 for sizes from L = 24 to 150. 
The curves collapse in one curve. 

Fig. |TT|: x K L~ lK / v " as function of U K with 7 K /z/ K = 1.74 for sizes from L = 24 to 150. 
The curves do not collapse in one curve. 

Fig. 12: rj K and r] M as function of critical temperature chosen. At] k and Ar] M are 



the estimate from a direct fit of \ as function of U. This give estimates of the critical 
temperatures in agreement with (p2|) and (^). Black circles show the result using dynamical 
properties. 

Fig. [14|: (k) as function of the time t for various system sizes at the critical temperature 
T C K = 0.5122. The curves for the biggest sizes L = 201 and L = 300 collapse within the 
errors. 



Fig. [13|: (k) as function of the time t for temperatures around T C K = 0.5122. A linear 



behavior appears only for T = T* 



Fig. [15]: Binder cumulant U K as function of time t at T C K = 0.5122. The slope gives d/z K . 
Only data between t = 300 and t = 10, 000 (shown by the arrows) is used for the fit. 

Fig. [16]: Binder cumulant U K and susceptibility \ K as function of V* at T C K = 0.5122. 
The slopes give v n and 7 K respectively. We use only data between t = 300 and t = 10, 000 
(shown by the arrows) for the fit. 
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Fig. [17]: \ K as function of the binder cumulant U K at T C K = 0.5122. The slope gives 
(2 — i] K )/d. We use only data between t = 300 and t = 10,000 (shown by the arrows) for 
the fit. 



Fig. [L8|: Phase diagram of the Ising-Xy model. Solid and dotted lines indicate continuous 
and first-order transitions respectively. The filled black circles show the possible fixed points. 
A and C are free parameter. Similar phase diagram appears for the J\ — J2 model [24|, the 



Zig-Zag model |TJ,|ID| and the RSOS model 

Fig. IS: Binder's parameter U M for the U(l) order parameter function of the temperature 
for various sizes L. The arrow shows the temperature of simulation T s = 0.511. The scales 



is similar to those of Fig 

„M 



Fig. ||: x M L 2 ~ v as function of U M with r/ M = 0.34 for the sizes L = 81, 105, 123 



and 150. The curves do not collapse in one curve. 



Fig. [2T 



X L v as function of U with rf 
and 150. The curves collapse in one curve. 

M L 2-v M as f unc tion of U M with 7] M 



Fig. |22 
and 150 



X 



0.375 for the sizes L = 81, 105, 123 



0.41 for the sizes L = 81, 105, 123 



The curves do not collapse in one curve 



Fig. 23: x K L 



as function of U K with 7 K /z/ K = 1.885 for the sizes from L = 24 to 



150. The curves do not collapse in one curve. 



Fig. H 
and 150 
Fig. 



X M L 2 ~ v '" as function of U M with r] 3 
The curves do not collapse in one curve. 

^: Helicity T for the U(l) symmetry as function of the temperature for various sizes 
L. The arrows show the temperature of simulation T s = 0.511, The critical temperatures 
for the U(l) symmetry T C A/ and the Ising symmetry T C K , To showing the best fit with the 
scaling form (0). The dashed line 2T/ir represents the universal ferromagnetic jump. 

Fig. ^: Helicity T for the U(l) symmetry as function of the lattice size L at T C M . The 
jump is smaller than the universal ferromagnetic jump (dashed line). 



0.461 for the sizes L — 81, 105, 123 



Fig. ^: x L 1 as function of U M with 77 = 0.22 for the sizes L = 81, 105, 123 
and 150. The curves do not collapse in one curve. 

Fig. |28|: (M) as function of the time t for various temperature. A linear behavior 
appears only for T = T C A/ = 0.5102. The deviation from a linear behavior is clearly visible 
for T = T = 0.5010. 
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